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ABSTRACT 

Context. The Gaia-ESO Public Spectroscopic Survey is obtaining high quality spectroscopy of some 100,000 Milky Way 
stars using the FLAMES spectrograph at the VLT. 

Aims. This paper describes the analysis of UVES and GIRAFFE spectra acquired in the fields of young clusters whose 
population includes pre-main sequence (PMS) stars. 

Methods. Both methods that have been extensively used in the past and new ones developed in the contest of the Gaia- 
ESO survey enterprise are available and used. The internal precision of these quantities is estimated by inter-comparing 
the results obtained by such different methods, while the accuracy is estimated by comparison with independent 
external data, like effective temperature and surface gravity derived from angular diameter measurements, on a sample 
of benchmarks stars. A validation procedure based on such comparisons is applied to discard spurious or doubtful 
results and produce recommended parameters. Specific strategies are implemented to deal with fast rotation, accretion 
signatures, chromospheric activity, and veiling. 

Results. The analysis carried out on spectra acquired in young clusters’ fields during the first 18 months of observations, 
up to June 2013, is presented in preparation of the first release of advanced data products. These include targets in 
the fields of the pOph, Ghal, NGC2264, 7 Vel, and NGG2547 clusters. Stellar parameters obtained with the higher 
resolution and larger wavelength coverage from UVES are reproduced with comparable accuracy and precision using 
the smaller wavelength range and lower resolution of the GIRAFFE setup adopted for young stars, which allows us 
to provide with confidence stellar parameters for the much larger GIRAFFE sample. Precisions are estimated to be « 

120K r.m.s. in T^a, ~0.3 dex r.m.s. in logg, and ~0.15 dex r.m.s. in [Fe/H], for both the UVES and GIRAFFE setups. 

Key words, open clusters and associations: general - surveys - methods: data analysis - stars: pre-main sequence - 
stars: fundamental parameters - open clusters and associations: individual: pOph, Chal, NGC2264, 7 Vel, NGC2547 


1. Introduction 

Spectrum analyses of pre-main sequence stars (PMS) re¬ 
quire special techniques, notably for dealing with peculiar¬ 
ities of cool, low-mass members of young clusters. Optical 
spectra of such stars may include the presence of veiling, 
large broadening due to fast rotation, emission lines due 
to accretion and/or chromospheric activity, and molecular 
bands. The subtraction of inhomogeneous and variable neb¬ 
ular emission may also be problematic and some residual 
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features can remain in spectra of some young clusters’ mem¬ 
bers after the sky-background removal. 


One of the main objectives of the Gaia-ESO Survey is to 
provide radial velocities (RV) with a precision « 0.2 - 0.25 
km s“® for stars in young open clusters, to complement Gaia 
proper motions wit h comparable accuracy for a statistically 
significant sample (Gilmore et al. 2012 Randich et al.|20l3 1 
reaching also fainter targets. 'This survey also complements 
Gaia by deriving metallicity and detailed abundances for 
several elements, including lithium, which is particularly 
relevant in the studies of the evolution of low-mass stars 
and in the determination of clusters’ age. This requires a 
derivation of all fundamental parameters (effective temper- 
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atures Teff, metallicity [Fe/H], surface gravity logg, and 
projected rotational velocity usini) independently of the 
Gaia results. 

The Ha profile of such young low-mass stars bears infor¬ 
mation on their chromospheric activity, accretion rate, and 
mass loss. Because of their common origin, strong accretion 
is expected to be correlated with veiling; this can be used 
for both checking our results, as no correlation would be 
indicative of large uncertainties, and exploring the extent 
and details of such a correlation. Chromospheric activity 
is known to depend on stellar rotation and both evolve in 
time; the Gaia-ESO is also going to provide the possibility 
of exploring the activity-rotation relation and their evolu¬ 
tion on a large sample of young stars. 

The Gaia-ESO target selection aims at producing un¬ 
biased catalogues of stars in open clusters. Selection cri¬ 
teria based mainly on photometry, supported, when possi¬ 
ble, by kinematic memberships, have been adopted for this 
purpose, although this implies that a large number of non¬ 
members are also observed, which are identified a posteriori 
from the radial velocity measurements (Bragaglia et al., in 
prep.). In our case, the GIRAFFE targets are late-type (F 
to early-M) stars in the magnitude range 12<V<19 mag, in 
the PMS or main sequence (MS) phase. Based on available 
information, the selection of UVES targets tries to include 
only slowly rotating {v sin i< 15 km s“^) single G-K stars in 
the magnitude range 9<V<15 without or with weak accre¬ 
tion (Macc< 10“^° Mq yr“^). To optimise the throughput 
of the survey, observations of cool stars in the fields of young 
open clusters are only carried out in the GIRAFFE/HR15N 
setup (R=17 000 , A from 6470 to 6790 A) and the Red 580 
UVES setup (R=47 000 centred at A = 5800 A with a spec¬ 
tral band of 2000A). The Medusa mode of the fibre fed 
system is used throughout the survey, allowing the simul¬ 
taneous allocation of 132 and 8 fibres feeding GIRAFFE 
and UVES, respectively, with about 20 (GIRAFFE) and 1 
(UVES) fibres used to observe the sky background spec¬ 
trum. 

The GIRAFFE/HR15N setup covers both Ha and 
Li (6707.84 A) lines, and it is therefore particularly use¬ 
ful for the study of young stars. However, T^g, log 5 , and 
[Fe/H] diagnostics in this wavelength range are poorer than 
in other settings and still not satisfactorily reproduced by 
theoretical models. For example, the paucity of Fe lines in 
the HR15N spectral range makes it difficult to derive both 
loggf and [Fe/H] in G-type stars from the analysis of the 
equivalent widths of Fei and Fell lines. 

This paper presents the analysis of the Gaia-ESO spec¬ 
tra in the fields of young open clusters (age < 100 Myr) 
and is one of a series presenting a description of the Gaia- 
ESO survey in preparation of its first release of advanced 
data products. The Gaia-ESO scientific goals, observations 
strategies, team organisation, target selection strategy, data 
release schedule, data reduction, analysis of OBA-type and 
FGK-type stars not in the fields of young open clusters, 
non-standard objects and outliers, external calibration, and 
the survey-wide homogenisation process are discussed in 
other papers of this series. 

The paper is organised as follows. In Sect.j^ the data 
analysed in the first two Gaia-ESO internal data releases 
are presented. In Sect.the principles and general strate¬ 
gies of the Gaia-ESO PMS analysis are outlined. Methods 
and validation for the initial raw measurements, fundamen- 
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tal parameters (Teff,logg, [Fe/H], micro-turbulence veloc¬ 
ity, veiling, and usini), and derived parameters (chromo¬ 
spheric activity, accretion rate, and elemental abundances) 
are presented in Sects. ilH and [^ The conclusions are in 
Sect. [3 


2. Data 

The survey’s analysis is performed in cycles, following the 
data reduction of newly observed spectra. Each new anal¬ 
ysis cycle improves upon the last one with updated input 
data (e.g. atomic and molecular data), improved analysis 
methods, and improved criteria to define the final recom¬ 
mended parameters. At the end of each cycle an internal 
data release (iDR) is produced and made available within 
the Gaia-ESO consortium for scientific validation. The de¬ 
scription of the methods and recommended parameters cri¬ 
teria given in this paper applies to the analysis of the first 
2 years of observations, which will form the basis of the 
first release of advanced data products to ESO. The vali¬ 
dation procedures presented in this paper consider the first 
18 months of observations (iDRl and iDR2). 

The young open clusters observed in the first 18 months 
of observations are listed in Table]^ along with the total 
number of observed stars for each cluster and the num¬ 
ber of stars identified as T Tauri from the properties of Ha 
emission, spectral type, and Li absorption (see Sect. 4.1). 
A total of 813 and 45 T Tauri stars have been identified in 
the GIRAFFE and UVES spectra, respectively. The mem¬ 
berships of these young clusters, including stars not clearly 
showing the T Tauri distinctive features, will be di scussed 
in other Ga ia-ESO science verification papers (e.g., Jeffries 
et al.|20i4 1. The 200 Myr cluster NG C6705 (Mil) has been 
observ ecTin different setups (see e.g. Gantat-Gaudin et al. 
2014b I to allow inter-comparison and validation of the anal¬ 


ysis methods across the survey and, for this purpose, is in¬ 
cluded in our analysis. 


3. General analysis strategy 

The Gaia-ESO consortium is structured in several work¬ 
ing groups (WGs). The analysis of PMS stars is carried 
out by the WG12, to which six nodes contributed: INAF- 
Osservatorio Astrofisico di Arcetri, Gentro de Astrofisica de 
Universidade do Porto (GAUP), Universita di Gatania and 
INAF-Osservatorio Astrofisico di Gatania (OAGT), INAF- 
Osservatorio Astronomico di Palermo (OAPA), Universi- 
dad Gomplutense de Madrid (UGM), and Eidgenossische 
Technische Hochschule Zurich (ETH). 

The main input to the Gaia-ESO PMS spectrum analy¬ 
sis consists of UVES and GIRAFFE spectra of cool stars in 
the field of young open clusters. The preliminary selection 
criteria are briefly outlined in Sect.[ 2 and will be detailed in 
one of the papers of this series (Bragaglia et al., in prep.). 
Th e data reduction is performed as described in [Sacco] 


et al. (2014) for the UVES spectra, and Lewis et al. (in 
prep.) for the GIRAFFE spectra. These are put on a wave¬ 
length scale and shifted to a barycentric reference frame. 
Sky-background subtraction, as well as a normalisation to 
the continuum, is also performed in the data reduction pro¬ 
cess. Multi-epoch spectra of the same source are combined 
in the co-added spectrum. Quality information is provided 
including variance spectra, S/N, non-usable pixels, etc. Ad¬ 
ditional inputs are the radial and rotational velocities, as 
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Table 1. Young open clusters (age < 100 Myr) observed by the Gaia-ESO survey in the first 18 months of observations, whose 
analysis is discussed in this paper. The cluster NGG6705 (Mil) has also been included for validation and comparison across the 
survey. 


Cluster 

Approximate Age 
(Myr) 

Distance 

(^Pc) 

All 

GIRAFFE 

WTTS/CTTS 

All 

UVES 

WTTS/CTTS 

iDR 

p Oph 

1 

120 

200 

30 

23 


5 

2 

Chal 

2 

160 

674 

93 

49 


14 

1 , 2 

NGC2264 

3 

760 

1706 

446 

118 


23 

2 

7 Vel 

~ 5-10 

350 

1242 

200 

80 


2 

1 , 2 

NGC2547 

35 

361 

450 

44 

26 


1 

2 

NGC6705 

250 

1877 

1028 

0 

49 


0 

2 


described in Gilmore et al. (in prep.) and Sacco et al. (20141, 
and photometric data. Clusters’ distance and reddening are 
also considered as input to the spectrum analysis validation. 

Double-lined spectroscopic binaries (SB2) and multiple 
systems are identified by looking at the shape of the cross¬ 
correlation function. These stars are excluded from the cur¬ 
rent analysis and a multiplicity flag is reported in the final 
database. 

To ensure the highest homogeneity as possible in the 
quantities derived, all the different Gaia-ESO spectrum 
analysis methods adopt the same atomic and molecular 
data (Heiter et al., in prep.), as well as the same s et of 
model atmospheres (MARCS, |Gustafsson et al.' 20081. 

The output parameters of the Gaia-ESO PMS spectrum 
analysis are listed in Table]^ To apply a detailed quality 
control on the output parameters and optimise the anal¬ 
ysis according to the star’s characteristics, these are di¬ 
vided into three groups: raw, fundamental, and derived. Raw 
parameters are the Ha emission and Li equivalent widths 
(VE(Ha) and W{h\)), and the Ha width at 10% of the line 
peak (Ha 10%, see, e.g., Natta et al. 20041. These are di¬ 
rectly measured on the input spectra and do not require 
any prior information. Their measurement is carried out 
before any other procedure to identify PMS stars and their 
values are used for optimising the evaluation of the funda¬ 
mental parameters in one of the methods used (see Sect.j^. 
Besides Teff, logg, and [Fe/HQ the fundamental parame¬ 
ters derived include also micro-turbulence velocity (^), pro¬ 
jected rotational velocity (usinz), veiling (r, see, e.g., |Har- 
tigan et al. 1988|), and a gr avity-sensitive spectral index 
( 7 , see Damiani et al. 2014|). Finally, the derived param¬ 
eters are those whose derivation requires prior knowledge 
of the fundamental parameters, i.e. elemental abundances 
(loge(X)Q, mass accretion rate (Mace), chromospheric ac¬ 
tivity indices (AH7(Ha)chr and AlP(H/3)chr), and chromo¬ 
spheric line fluxes (F(Ha)chr and E(H/3)chr)- 

Most parameters listed in Table|^are derived from both 
UVES and GIRAFFE spectra, with the exception of 
AkP(H/3)chr, P’(H/3)chr, and loge(X), that are derived from 
UVES spectra only, and the gravity-sensitive spectral index 
7 , which is derived from the GIRAFFE spectra only (see 
Sect. I . 

In general, whenever possible, the same parameter is de¬ 
rived by different methods; this allows a thorough check of 
the derived parameters by inter-comparing the results and 


^ The solar Fe abundance of Grevesse et al. (20071, log e(Fe )0 = 
7.45, is adopted. 

^ log e(X) = log[A(X)/A(H)] -f 12, i.e. a logarithmic abundance 
by number on a scale where the number of hydrogen atoms is 
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flagging discrepant results, which are then used to outline 
possible weaknesses of the methods and discard unreliable 
results. In the absence of significant biases, the results from 
different methods are combined taking a cr-clipped aver¬ 
age to obtain the recommended parameters. If significant 
biases are present, all results obtained with a method that 
can give rise to inaccurate or unreliable results in some 
ranges of parameters are rejected before combining the re¬ 
sults as above. These general criteria, whose application 
is discussed in details in Sects. |4}|^ are firstly applied on 
the raw parameters, then on the fundamental parameters, 
and Anally on the derived parameters. The fundamental pa¬ 
rameters are also validated by comparing the results of the 
analysis methods applied to our spectra against fundamen¬ 
tal parameters from angular diameter and parallax mea¬ 
surements (Sect. 5 . 61 . A comparison with derived from 
photometry for objects that are not affected by photometric 
excesses is reported in Appendix[C| The recommended raw 
and fundamental parameters are then used to produce the 
recommended derived parameters. When satisfactory com¬ 
parisons cannot be achieved, recommended parameters are 
not provided and only results from individual nodes are 
made available. Recommended parameter uncertainties are 
estimated as both node-to-node dispersion and as average of 
individual node’s uncertainties. The final results minimise 
- as much as possible - biases that can affect individual 
methods and the associated uncertainties take differences 
that may arise from the use of different methods and algo¬ 
rithms into account. 

Final results are further validated by a general analysis 
of the output log g-Tgg diagram, consistency of the param¬ 
eters, and overview of the results based on the comparison 
of different clusters. 


4. Raw measurements 

Measuring the raw parameters before carrying out any 
other analysis allows us to: (a) identify stars with strong 
accretion whose spectra may be affected by veiling; (b) per¬ 
form a quality control on the raw parameters before they 
are used in the subsequent analysis; and (c) apply the ap¬ 
propriate masks to the spectra for the determination of fun¬ 
damental parameters. 

To derive raw parameters from a large dataset of spectra 
it is convenient to use procedures that are as automatic as 
possible. However, in the case of PMS sources extending to 
M spectral type, such procedures must also be capable of 
dealing with large rotational broadening and the presence 
of molecular bands. Here different methods are used, with 
different levels of automatism, which allows to examine the 
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Table 2. Gaia-ESO PMS analysis output parameters. Columns 2-13 list the number of stars in each cluster for which the parameter 
was derived from GIRAFFE (G) and UVES (U) spectra separately in iDR2. Lithium parameters and usini counts include upper 
limit estimates. Accretion and chromospheric activity parameter counts include only non negligible values. For the elemental 
abundances, the maximum number of derived values for each star/element is reported. See text for an explanation of the notation 
used. 


Parameter 

p Oph 

Chal 

NGC2264 

7 Vel 

NGC2547 

NGC6705 


G 

U 

G 

U 

G 

U 

G 

U 

G 

U 

G 

U 

raw 

lU(Ha) 

25 

5 

87 

14 

387 

24 

203 

2 

106 

1 

0 

0 

lU(Li) 

189 

23 

633 

47 

1610 

114 

1186 

75 

404 

25 

708 

48 

Ha 10% 

33 

5 

103 

14 

807 

23 

264 

2 

239 

1 

0 

0 

fundamental 

TeS 

170 

21 

572 

39 

1324 

70 

1104 

51 

361 

24 

394 

32 

logff 

170 

21 

156 

39 

226 

70 

350 

51 

106 

24 

150 

32 

7 

156 


508 


1199 


1043 . 


337 


0 


[Fe/H] 

170 

21 

515 

39 

1203 

70 

1018 

51 

311 

24 

360 

32 



14 


23 


42 


46 


15 


30 

tjsinz 

154 

23 

521 

42 

1192 

83 

1004 

75 

332 

25 

107 

33 

r 

4 

3 

20 

7 

77 

6 

5 

0 

5 

0 

0 

0 

derived 

loge(Li) 

154 

23 

514 

40 

1203 

80 

1017 

57 

311 

25 

356 

31 

loge(A) 


15 


28 


39 


46 


14 


31 

Macc 

14 

4 

56 

7 

212 

11 

40 

1 

21 

0 

0 

0 

AlF(Ha)ehr 

21 

12 

69 

29 

267 

50 

205 

18 

115 

16 

61 

0 

AlU(H/3)ehr 


10 


18 


42 


14 


12 


0 

F(Ha)chr 

21 

12 

65 

28 

265 

47 

199 

17 

105 

16 

47 

0 

E(H/3)chr 


10 


17 


41 


14 


12 


0 


presence of biases, eliminate systematic discrepancies, and 
combine the results with a tr-clipping to disregard casual 
mistakes and outliers. 

In the following we briefly describe the methods used to 
derive the raw parameters. 


4.1. Ha equivalent width and Ha width at 10% of the line 
peak 


Spectra with Ha in emission are examined to identify stars 
with strong accretion, and therefore likely to be affected by 
veiling, using their W (Ha) and Ha 10% measurements. 

The Arcetri node measures W (Ha) and Ha 10% on the 
continuum-normalised co-added spectra of all stars that 
clearly show Ha emission, using a semi-automatic proce¬ 
dure. After manually defining the wavelength range and 
level of continuum, IT(Ha) is calculated by a direct inte¬ 
gration of the flux above the continuum, while Ha 10% is 
derived by considering the level corresponding to 10 % of the 
maximum flux above the continuum in the selected wave¬ 
length range. All measurements are visually checked and 
repeated in case of miscalculation (e.g. due to the presence 
of multiple peaks). Uncertainties are estimated using multi¬ 
epoch observations of stars belonging to the first two young 
clusters that have been observed (i.e. 7 Vel and Chal). 
Specifically, W (Ha) and Ha 10% are first measured on each 
spectrum before co-adding, then the relative uncertainty for 
each star is estimated as AW = 2 \ Wi — H 2 I /(lUi + W 2 ), 
where Wi and W 2 are two measurements for the same star 
from spectra observed at different epochs. A similar formula 
is used for AHa 10%. Finally, the median of AW (Ha) and 
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AHa 10% are assumed as the relative uncertainties for all 
star^ 

The CAUP node makes use of an automatic idl| 3 pro¬ 
cedure to first select stars with Ha in emission and then 
measure Ha 10% and W (Ha) on the normalised spectra. 
Measurement uncertainty is evaluated from the spectrum 
S/N. 

The OACT node pre-selects spectra with Ha in emis¬ 
sion by visual inspection. Then, W (Ha) and Ha 10% are 
measured using an IDL procedure. IT(Ha) is measured by 
direct integration of the Ha emission profile and its uncer¬ 
tainty evaluated by multiplying the integration range by 
the mean error in two spectra regions close to the Ha line. 
Ha 10% uncertainty is evaluated by assuming an error of 
10 % in the position of the continuum level. 

The PA PA node employed two me thods, one based on 
DAOSPEC (Stetson & Pancino 20081 and an IDL proce¬ 
dure, the other on a combination of IRAF and IDL tools. 
In the first method, DAOSPEC is used to perform a con¬ 
tinuum fit of the spectral region around Ha. The Ha profile 
in the unnormalised input spectrum is masked by giving as 
input to DAOSPEC a variable FWHM that takes the ro¬ 
tational and instrumental profiles into account. The fitted 
continuum is then used to normalise the input spectrum. 
Such a continuum normalised spectrum is used to measure 
Ha 10% by an automatic IDL procedure. Since the uncer¬ 
tainties are assumed to be dominated by the fitting of the 
continuum, this is repeated four times using different orders 
(10, 15, 20 and 25) of the polynomial fitting in DAOSPEC. 
The resulting Ha 10% values are then averaged to produce 
the final result. In the second method the normalisation is 


^ This may also be linked to Ha variability. 

^ IDL(r) (Interactive Data Language) is a registered trademark 
of Exelis Visual Information Solutions. 
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performed through IRAF with three different orders of the 
polynomial fitting (2, 5, and 10), then W(Ha) and Ha 10% 
are measured with an automatic IDL routine and uncer¬ 
tainties derived as above. A final visual inspection is per¬ 
formed to check the results and identify broad emission and 
P Cygni-like profiles. For the first data release both meth¬ 
ods where used, while in the second data release only the 
second method was used. 


The final spectra of NGC2264 are affected by some 
residual nebular emission, and a good subtraction of this 
contribution to the Ha emission line cannot be achieved as 
the nebular emission is concentrated in the region near the 
Ha line peak and is spatially variable (see a detailed de¬ 
scription of this topic for the a nalogous case of the cluster 
NGC6611 in Bonito et al.|2013 l. In this case, additional vi¬ 
sual inspection of the spectra was necessary to ensure that 
the narrow nebular emission does not affect significantly 
the measurements. 


• s ’i* 

• X • • .* . 

• ... •. » *» . 



400 

Ho 10% (km/s) 


p Oph (1 Myr) O • 
Cha I (2 Myr) 
NGC2264 (3 Myr) 0« 
y Vel (5-10 Myr) 
NGC2547 (35 Myr) Ai 
. . I 

600 


Fig. 1. IF(Ha) vs. Ha 10% for all young open clusters observed 
in the first 18 months of observations. Filled symbols are used 
for stars classified as CTTS, open symbols for stars classified as 
WTTS. 


In the node-to-node comparison of the VF(Ha) results, 
average differences and dispersions ^ 5 A were found in the 
analysis of both UVES and GIRAFFE spectra, with only a 
few outliers. Average differences in Ha 10% in the node-to- 
node comparison was ^10 kms“^, with a dispersion ~50 
kms“^. 


Only a la-clipping was therefore applied before com¬ 
puting the average kF(Ha) and Ha 10% as recommended 
values. The recommended uncertainty was given, conser¬ 
vatively, as the largest amongst the average of individual 
uncertainties and the standard deviation of the mean. 


The recommended Ha 10% is used, together with the 
recommended W(Li) (Sect. 4.21, in our WTTS/GTTS clas¬ 
sification. If the Ha is in emission and H^(Li)>100 mA, 
the star is identified as a TTauri. Following [White iS] 


Basri (20031, the TTauri star is then classified as GTTS 
if Ha 10%>270 kms“^. 

A comparison of kF(Ha) vs. Ha 10% for all young open 
clusters observed in the first 18 months of observations is 
shown in Fig.[l] Note that the correlation of the two param¬ 


eters is as expected from other works (e.g.. White & Basri 


20031 and the fraction of GTTS consistently decreases with 
the age of the cluster. 


4.2. Li equivalent width 


At young ages, the Li doublet line is often in the saturation 
regime and the rotational broadening frequently dominates. 
As a consequence, in general, a direct profile integration 
of the Li line is to be preferred to a Gaussian or a Voigt 
profile fitting in deriving WilA). Furthermore, due to ro¬ 
tational broadening, the integration wavelength interval is 
very different from one spectrum to another. The Li dou¬ 
blet is also superimposed to molecular bands in spectra of 
M-type stars, which makes the placement of the continuum 
difficult, particularly when using automatic procedures. In 
such cases, despite being slow and prone to human error 
and subjective choices, interactive procedures, like those 
available in IRAF, remain one of the best options for mea¬ 
suring VF(Li), at least for comparison purposes. Weak Li 
lines in slow-rotating stars, on the other hand, can reliably 
be fitted with a Gaussian or a Voigt profile and integrated 
analytically, a method that can be easily implemented in 
automatic procedures and is more accurate than low-order 
numerical integration at low S/N. 

The Gaia-ESO PMS analysis makes use of three inde¬ 
pendent methods for deriving W{h\) from the GIRAFFE 
spectra: the direct profile integration available in the 


IRAF-splot proc edure (OAGT node); DAOSPEG (Stetson 
& Pancino 2008 OAPA node); and a semi-automatic IDL 
procedure specifically developed for the Gaia-ESO by the 
Arcetri node. 

The IRAF-splot task was applied by the OAGT node 
to the unnormalised spectra to make use of the built-in 
Poisson statistics model of the data. Such measurements 
are only performed in those cases where the Li line and 
the nearby continuum are clearly identifiable, which implies 
that, in general, small IF(Li) (< 10 niA), low S/N (< 20), 
and spectra with very high usinf (> 200 kms“^) are not 
considered. 

The DAOSPEG (OAPA) measurements were applied to 
all iDRl spectra with S/N> 20, and spectra with S/N< 20 
showing a strong lithium line. The spectra have been re¬ 
normalised prior to the equivalent width determination us¬ 
ing high order Legendre polynomial fitting, which allows to 
follow the shape of molecular bands in M-type stars still 
maintaining a good agreement with the continuum of ear¬ 
lier type stars. The typical width of absorption lines in 
each spectrum has been estimated by convolving the in¬ 
strumental and rotational profile using v sini from the data 
reduction pipeline. Relative internal uncertainties are al¬ 
ways better than 5% for large equivalent widths (> 200 mA) 
and degrade up to ^ 50% for very low equivalent widths 

(^ 10 mA). 

The semi-automatic IDL procedure developed by the 
Arcetri node performs a spline fitting of the continuum over 
a region of ±20 A around the Li line using an iterative a- 
clipping, and masking both the Li line and the nearby Gai 
line at 6717.7 A. When the automatic continuum fitting is 
not satisfactory (generally for poor S/N spectra or M-type 
stars), the fit is repeated by setting manually the continuum 
level. The W (Li) is then computed by direct integration of 
the line within a given interval, which depends on the stellar 
rotation and was determined by measuring the line widths 
on a series of rotationally broadened syn thetic spectra. Er¬ 
rors are derived using the Gayrel (1988) formula; when no 
Li line (including blends) is visible, the upper limit is set 
as three times the error. 
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Fig. 2. Illustrative node-to-node W{\A) comparison for 7 Vel. Left panel: Comparison between OAPA (DAOSPEC, iDRl) and 
Arcetri (iDR2) for GIRAFFE spectra. Right panel: Comparison between Arcetri code and OACT (IRAF) for UVES spectra. 
Arrows indicate upper limits. 


The contribution of lines blended with Li in the GI¬ 
RAFFE spectra was estimated, after the determination of 
the fundamental parameters (Sect.l^, by a sp ectral syn¬ 
thesis using Sp ectroscopy Made Easy (SME, Valenti & 
Piskunov|1996 1 with MARCS model atmospheres as input, 
taking the star’s logf/, and [Fe/H] into account. For 
solar metallicity dwarfs a bove 4000 K the estimat ed blends 
are in agreement with the Soderblom et al. ( 1993[ | relation. 

Four nodes (Arcetri, CAUP, OACT, and UCM) calcu¬ 
lated W (Li) in the UVES spectra. At the UVES resolution, 
when the star is slow rotating (usinf< 25 kms“^) and the 
S/N is sufficiently high (S/N> 60), it is possible to de-blend 
the Li line from the nearby features. 

Both the CAUP and UCM nodes employed the splot 
task in IRAF on the unnormalised UVES spectra. When 
the Li line and the nearby blends, mainly with Fe lines, 
are distinguishable, these are de-blended, in which case a 
Gaussian fitting to the line profile is adopted. On the con¬ 
trary, when the lines are indistinguishable, the blends con- 
tribut ion is estimat ed using the ewfind driver within MOQG 
code (Sneden 19731, and a direct integration of the line is 
adopted! 

The Arcetri node adopted the same method used for 
GIRAFFE (see above), except in those cases where it was 
possible to de-blend the line using IRAF as done by the 
CAUP and UCM nodes. When this was not possible, the 
blends were estimated using SME. 

For iDRl, the OACT node employed IRAF as for the 
GIRAFFE spectra, using SME to estimate the blends. For 
iDR2, IV(Li) was derived by spectra subtraction with the 
template having the closest fundamental parameters but no 
(or negligible) Li absorption. In this latter case the blends 
are removed by the spectra subtraction itself. 

It is worth stressing that the PMS analysis output 
includes both blends-corrected and -uncorrected tV(Li). 
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When a node does not provide blends-corrected lU(Li), 
this is estimated using SME and the node’s fundamental 
parameters if available. Note that, in the analysis of GI¬ 
RAFFE spectra, blends are estimated using SME in all 
cases; the recommended blends-corrected W (Li) are calcu¬ 
lated from the recommended blends-uncorrected W (Li) us¬ 
ing the recommended fundamental parameters. Conversely, 
in the analysis of UVES spectra the recommended blends- 
corrected lU(Li) are derived by averaging the node values, 
as discussed below. 

The blends-uncorrected VF(Li) obtained by the three 
different methods from GIRAFFE spectra were first com¬ 
pared to check for systematic differences before combining 
them to produce the final results (see Fig. II for an illustra¬ 
tive comparison). After discarding results of one node not 
consistent with the other two, no significant bias remained 
and the relative standard deviation of the W difference was 
at 20% level. Also, no trend of the node-to-node differences 
with respect to S/N nor usinf was present in the selected 
results. As a conservative uncertainty estimate on the rec¬ 
ommended W (Li) we adopted the larger of the standard 
deviation and the mean of the individual method uncer¬ 
tainties. In iDRl, 90% of the IV(Li) measurements have 
relative uncertainties better than 14% and 28% in 7 Vel and 
Chal, respectively, the differences being due to the higher 
fraction of stars of low T^g and spectra with lower S/N in 
Chal with respect to 7 Vel. About 90% of all the iDR2- 
GIRAFFE IV(Li) measurements have uncertainty better 
than 16 mA. 

In the W (Li) UVES measurements no systematic devi¬ 
ation nor trends of the node-to-node differences with S/N 
nor V sin i were found from the node-to-node comparison 
(see Fig.H for an illustrative comparison) and the recom¬ 
mended values were derived by taking the mean with a Icr- 
clipping. In iDRl the median uncertainties are 3 mA (4%) 
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and 10 mA (3%) for 7 Vel and Chal, respectiveljjj About 
90% of all the iDR2-UVES WilA) measurements nave un¬ 
certainty better than 22 mA. 

When all W{h\) measurements for a given star are 
flagged as upper limits, the recommended W{L\) is also 
flagged as an upper limit and the lowest measurement is 
adopted. Conversely, when at least one W{Li) measure¬ 
ment for a given star is not flagged as upper limit, all upper 
limit estimates for that star are disregarded and the recom¬ 
mended value is derived as above. 



Fig. 3. Blends-corrected W{L\) vs. Tee for all young open clus¬ 
ters observed in the first 18 months of observations. For com¬ 
parison, the lower and upper envelope in the Pleiades are shown 
as solid lines, while dashes and dotted lines are the upper en¬ 
velopes of IC2602 and the Hyades, respectively. Dots are used 
for non-members. Upper limits are not included for clarity. 


Recommended blends-corrected lU(Li) vs. are shown 
in Fig.[^ compared with the Pleiades upper and lower en¬ 
velope a nd the upper envelopes of IC2602 (30 Myr) and the 
Hyades (|Sestito &: Randich| 2005 and references therein). 
Note the increasing Li depletion with age at Teff> 3500 K 
and the lack of depletion at lower Teff as expected. A com- 
parison with theoretical models can be found in | Jeffries] 


et al. (20141 for 7 Vel. The NGC2547 Li depletion pattern 
is found in remarkable agreement with jJeffries & Oliveira] 
( 200 ^. 


5. Fundamental parameters 


Two nodes (OACT and OAPA) provide fundamental pa¬ 
rameters from the analysis of GIRAFFE spectra, and four 
nodes (Arcetri, CAUP, OACT, and UCM) from the analysis 
of UVES spectra. With the exception of OACT, the nodes 
analysing UVES spectra use similar, classical pro cedures, 
i.e. measurement of equivalent widths and MOOG (Sneden 


19731, enforcing the usual equilibrium relations. However, 
different strategies are adopted for the selection of the lines 
to be used and in the automatisation of the procedures, as 
described in Sects. (53115.51 

As anticipated in Sect.|^ the validation of fundamental 
parameters is carr ied out internally by a node-to-node com¬ 
parison (Sect. 5.71, and externally by comparisons with Teff 


® Note that, given the small number of measurements In UVES 
compared to GIRAFFE, the median uncertainty is chosen to 
characterise the internal precision achieved rather than the dis¬ 
tribution. 


and logg derived from angular diameter me asur ements on 
a sample of stars taken as benchmark (Sect. 5.6 1 


5.1. OACT 


A code that has been extensively used for determini ng fun¬ 
damental p arameters in PMS stars is ROTFIT (e.g., Frasca 
et al.| 200 ^ , which compares the target spectrum with a set 
of template spectra from ELODIE observations of s lowly- 
rotating non-active stars ( Prugniel fc Soubiran|200T I artifi¬ 
cially rotation-broadened and veiled at varying v sin i and r. 
In the following we report a brief summary of the method 
implemented by ROTFIT, together with some adaptations to 
the case at hand. 

In ROTFIT, the template spectra that most closely re¬ 
produce the target spectrum when broadened and veiled 
are selected, and their weighted average logg, and 
[Fe/H] assigned to the target star. As a figure of merit the 

calculated on the target spectrum and the rotational- 
broadened and veiled template spectrum is adopted. The 
weight used in the average is proportional to A discus¬ 
sion on the ROTFIT templates, together with the homogeni¬ 
sation with the Gaia-ESO spectrum analysis, is reported in 
AppendixfAl The ROTFIT analysis requires different wave¬ 
length masks for different type of objects. The masking cri¬ 
teria for the application to the Gaia-ESO survey are re¬ 
ported in AppendixlB] 

For the GIRAFFE/HR15N spectra the whole spectral 
range from 6445 to 6680 A and the ten best templates (i.e. 
with the lowest x^) ^re considered. The analysis of the 
UVES 580 spectra is independently performed on segments 
of 100 A each (excluding the parts containing strong tel¬ 
luric lines and the core of Balmer lines), by considering 
only the best five templates for each segment. The final 
parameters Teff, logg, [Fe/H], and usini are obtained by 
taking the weighted averages of the mean values for each 
segment, with the weight being proportional to x~^ and to 
the amount of information contained in the segment, quan¬ 
tified by the total line absorption = ^{F\/Fq — l)(iA 
(where F\ / Fc is the continuum normalised flux). 

The Teff, log g, [Fe/H], and v sin i uncertainties are given 
as the standard errors of the weighted means, to which the 
average uncertainties of the templates’ stellar parameters 
are added quadratically. These are estimated to be ±50K, 
iO.ldex, iO.ldex, and ±0.5kms“^ for Teff, logf;, [Fe/H], 
and usini, respectively. The target’s spectral type corre¬ 
sponds to that of the best template. 

The code also provides an estimate of the veiling by 
searching for the r value that gives the lowest x^. The de¬ 
termination of the fundamental parameters for a star with 
veiling is more uncertain than in the non-veiled case be¬ 
cause: (a) lines’ and molecular bands’ depth are smaller in 
veiled spectra; and (b) the determination of the veiling pa¬ 
rameter r implies the introduction of an additional degree 
of freedom in the parameters’ fitting, degrading the over¬ 
all accuracy with respect to the non-veiled stars. However, 
a preliminary identification of stars with mass accretion, 
whose spectra are expected to be affected by veili ng, can 
be done based on the values of Ha 10% or W{Cia) (White 


& Basri 2003 Natta et al. 


20041. It is therefore possible 


to restrict the veiling calculation to likely accreting stars 
only, thus preserving the ac curacy for stars for w hich no 
veiling is expected. Following White & Basri (20031, we as- 
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sume that stars with Ha 10%>270 kms“^ can be optically 
veiled, with an extr a margin to take uncertainties into ac¬ 
count (see Sect. 4.1 1 . 

Within the Gaia-ESO analysis, vsini is also provided 
by WG 8 for GIRAFFE spectra (Koposov et ah, in prep.), 
together with radial velocities and a first estimate of funda¬ 
mental parameters, using a conceptually similar approach 
but with a different fitting strategy and different templates. 
The comparison between the results of these two methods 
turned out to be useful in identifying WG 8 unsuccessful fit¬ 
ting for some stars with strong veiling and emission lines. 
An illustrative comparison o f the results of thes e two meth¬ 
ods for 7 Vel can be found in Frasca et al. (20141, who found 
a mean difference of « 2.88 km s“^ and ct « 6.27 for stars 
in the field of 7 Vel. An investigation on the lower limit 
imposed by the resolution of the instruments by means of 
Monte -Garlo simulations is also reported in [Frasca et al.| 
(20141, according to which the lower usini limit has been 
set to 7 kms“^ in GIRAFFE spectra and 3 kms“^ in UVES 
spectra. 


5.2. OAPA 

An alternati ve approach for GIRA FFE/HR15N spectra, 
proposed by Damiani et al. (20141, is based on spectral 
indices in different wavelength ranges of the spectrum. The 
derived spectral indices are calibrated against stars with 
known parameters, yielding quantitative estimates of T^e, 
logg, and [Fe/H]. This type of approach is usually applied 
to spectra with lower resolution than GIRAFFE. These 
narrow-band indices are affected by fast rotation: Teff be¬ 
comes unreliable for usinf>90 kms“^, [Fe/H] above 70 
kms“^, and log^ above 30 kms“^. Therefore, depending 
on the V sin i of the star, not all parameters can be pro¬ 
vided. Using an appropriate combination of flux ratios, this 


for details). 


5.3. Arcetri 

The Arcetri nod e adopted DOOp (DAO SPEG Option Op- 
timiser pipeline, Gantat-Gaudin et al. 2014a) for measur- 
ing line equivalent widths and FAMA (Magrini et al. 2013) 
for determining the fundamental parameters. Line equiva- 
lent widths are measured using Gaussian fitting after a re¬ 
normalisation of the continuum; W in the range between 
20-120 mA, for the Fei and Fell lines, and in the range 
between 5-120 mA, for the other elements, were used. The 
FAMA code makes use of the Fe i and Fe ii equivalent widths 
to derive the fundamental stellar parameters. The stellar 
parameters are obtained by searching iteratively for the 
three equilibria (excitation, ionisation, and trend between 
logn(Fei) and log(VF/A)), i.e. with a series of recursive 
steps starting from a set of initial atmospheric parameters 
and arriving at a final set of atmospheric parameters which 
fulfils the three equilibrium conditions. The convergence cri¬ 
terion is set using the information on the quality of the W 
measurements, i.e. the minimum reachable slopes are linked 
to the quality of the spectra, as expressed by the dispersion 
around the average < logn(Fei)>. This is correct in the 
approximation that the main contribution to the disper¬ 
sion is due to the error in the W measurement rather than 


to inaccuracy in atomic parameters, as e.g. the oscillator 
strengths (log 5 /). 


5.4. CAUP 

The fundamental parameters are automaticall y determined 
with a met hod used in previous works (e.g. Sousa et al. 

I now adapted to the Gaia-ESO survey. The 
jased on the excitation and ionisation balance 


2008 2011 


method is 

of iron lines using [Fe/H] as a proxy for the metallicity. 
The iron lines constraining the parameters were selected 
fr om the Gaia-ESO l ine list using a new procedure detailed 
Sousa et al. (2014). The equivalent widths are autom ati- 
ly measured using t he ARE^code (Sousa et al.|2007 1 fol- 


ca 


lowing the approach of Sousa et al. (2008 201 ip that takes 
the S/N of each spectrum into account. The stellar param¬ 
eters are computed assuming LTE using the 2002 version of 


MOOG (Sneden 


models. For t 


1973p and the MARGS grid of atmospheric 
the interpolation code provided 


tiis purpose 

with the MARGS grid was modified to produce an output 
model readable by MOOG. Moreover, a wrapper program was 
implemented to the interpolation code to automatise the 
method. The atmospheric parameters are inferred from the 
previously selected Fei-Feii line list. The downhill simplex 
(Press et al. 1992) minimisation algorithm is used to find 
the best parameters. In order to identify outliers caused by 
incorrect W values, a 3cr-clipping of the Fe I and Fe II lines is 
performed after a first determination of the stellar parame¬ 
ters. After this clipping, the procedure is repeated without 
the rejected lines. The uncertainties in t he stellar parame- 
ters are determined as in previous works ([Sousa et al.|2008 


5.5. UCM 


The UGM node employed the StePar code (Tabernero et al. 


2014 


2012 

2013 


2013), modified to operate with the spherical and 


non-spherical MARGS models. For iDRl the W measure¬ 


ment s were carried out with the ARES code (Sousa et al. 
2007|P For iDR2 the UGM node adopted the TAME code 


[Kan^ & Lee 2012 ® The manual mode has an interface 
that allows the user control over the W measurements to 
verify problematic spectra when needed. The [Kang fc Lee| 
(2012) approach to adjust the TAME continuum cr rejection 


parameter according to the S/N of each spectrum was fol¬ 
lowed. The StePar cod e computes th e stellar atmospheric 
parameters using MOOG (Sneden 1973). The 2002 and 2013 
versions of MOOG were used in iURl and iDR2, respectively. 
Five line lists were built-up for different regimes: metal rich 
dwarfs, metal poor dwarfs, metal rich giants, metal poor gi¬ 
ants and extremely metal poor stars. The code iterates un¬ 
til it reaches the excitation and ionisation equilibrium and 
minimises trends of abundance vs. log( IF/A). The down¬ 
hill simplex method (Press et al. 1992|) was employed to 
minimise a quadratic form composed of the excitation and 
ionisation equilibrium conditions. The code performs a new 
simplex optimisation until the metallicity of the model and 
the iron abundance are the same. Uncertainties for the stel- 


° ARES is available for download at 
http: / / WWW. astro.u p.pt/~sousasaK/are s/ 

' The approach of|Sousa et al. (2008) to adjust the ARES pa¬ 
rameters according to the S/JN of each spectrum was followed. 

® TAME is a tool that can be run in automated or manual mode. 
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lar paramet ers are derived as described in |Tabernero et al^ 
(2012 20131. In addition, a 3cr rejection of the Fe i and Fe ii 
lines is performed after a first determination of the stel¬ 
lar parameters; StePar is then re-run without the rejected 
lines. 


5.6. Comparison with benchmark stars 


The precision of the fundamental parameters can be as¬ 
sessed by comparison with results from accurate indepen¬ 
dent methods like interferometric angular diameter mea¬ 


surements (e.g. Boyajian et al. 2012a|b I which, in combi¬ 
nation with the Hipparcos parallax and measurements of 
the star’s bolometric flux, allow the computation of abso¬ 
lute luminosities, linear radii, and effective temperatures. 
As part of the Gaia-ESO activities, and also in support of 
the Gaia mission, a list of stars with accurate fundamen¬ 
tal parameters d erived from such independent methods is 
being compiled (Jofre et al. 2014 Heiter et al., in prep.) 
and included in the Gaia-ESO target list. For the range of 
parameters of interest to the PMS analysis, however, only 
very few benchmark stars spectra are available in iDRl and 
iDR2. 

A comparison of the iDR2-GIRAFFE fundamental pa¬ 
rameters of benchmark stars with those compiled from the 
literature is shown in Fig.|^ and Table[^ in the range of 
interest. In this case the T^s deviations are mostly within 
« 200 K. There is a systematic large deviation of OAGT val¬ 
ues above 6000 K. At lower temperatures, deviations larger 
than « 200 K are found in OAPA results for HD 10700. 
Therefore, although the sample analysed is limited, good 
results are found for both nodes, except for OAGT above 
6000 K. Excluding the OAGT values above this limit, the 
standard deviation is « 120 K for both datasets. The OAGT 
Teff upper limit for the GIRAFFE analysis was further low¬ 
ered to 5500 K based on compari son with Teff from photom¬ 
etry (Appendix!^ see also Sect. 5.7l. 

Deviations as large as almost 0.7 dex in log(; are found 
in the comparison with the benchmarks, with standard de¬ 
viation « 0.3 dex for both datasets. From the comparison 
with benchmarks alone it is not possible to identify a range 
in which one method performs better than the other. Indeed 
the node-to-node comparisons for each cluster outlined a 
rather complex situation that leads to the parameters se¬ 
lection described in Sect. 15.71 

In the parameters range of interest (i.e. excluding very 
metal-poor stars), [Fe/H] is approximately reproduced with 
a maximum deviation of 0.3 dex and a standard deviation 
of « 0.15 dex. 


The comparison of the UVES fundamental parameters 
of benchmarks with those compiled from the literature is 
shown in Fig.[^ and Table[3| The results for the solar spec¬ 
trum are outlined in TableU 

In general Teff deviations from benchmarks are all 
within 300 K (maximum) with a few outliers. Amongst 
these, UGM T^g for 61GygA differs by about 800 K, but 
this large deviation does not point to particular problems 
in some parameters’ range as verified through the node-to- 
node comparison. 

The OAGT systematic discrepancies in Teff above 
6500 K, on the other hand, indicate a T^g upper limit also 
for the validity of ROTFIT UVES analysis. This discrepancy 
is seen also in the node-to-node comparison, on which we 


estimate an upper limit of 6200 K for the validity of the 
OAGT results. 


Excluding such outliers, the standard deviation is « 
100 K for the Arcetri, GAUP, and UGM nodes, and « 120 K 
for the OAGT node, with average difference of 34 K for 
GAUP, 38 K for OAGT, and 55 K for Arcetri and UGM. 

Within the UVES dataset of young open clusters, very 
few sources have T^g <4000 K. In this range recommended 
data are based on OAGT results only, as the presence of 
molecular bands prevents to carry out analysis based on 
MOOG. 


The agreement in logg is approximately at the same 
level for all nodes. Benchmarks’ logg is generally repro¬ 
duced within a maximum deviation of « 0.7 dex and a 
standard deviation of « 0.3, only one Arcetri value deviat¬ 
ing more than that. 

[Fe/H] is generally reproduced within a maximum devi¬ 
ation of « 0.3 dex , except one and two measurements by 
the Arcetri and OAGT nodes, respectively, with deviations 
of « 0.5 dex. The standard deviation is < 0.1 dex for the 
GAUP and UGM nodes, « 0.2 dex for the Arcetri node, 
and « 0.3 for the OAGT node. The OAGT node tends to 
overestimate (underestimate) the metallicity below (above) 
[Fe/H]=0. However, this does not lead to significant system¬ 
atic differences in individual clusters and the OAGT results 
are therefore maintained. The node-to-node comparisons 
for individual clusters show that the Arcetri node systemat¬ 
ically overestimates [Fe/H[, which is not evident in the com¬ 
parison with the benchmarks possibly because of the large 
and coarse parameters’ distribution of this latter. To over¬ 
come this systematic behaviour, in iDR2 the value obtained 
by the Arcetri node for the solar spectrum ([Fe/H[=0.09, 
see Table|^ is subtracted in all measurements before com¬ 
puting the recommended [Fe/H[. The recommended [Fe/H[ 
agrees with the benchmarks within 0.15 dex r.m.s. 

Solutions with large uncertainties or large C 2 
kms“^) are disregarded by the nodes. Differences in ^ with 
respect to the values tabulated for the benchmarks are gen¬ 
erally below 1 kms“^. 

The recommended fundamental parameters are there¬ 
fore computed taking the average of the nodes’ results with 
a Itj-clipping when at least 3 values are provided. As dis¬ 
cussed above, below 4000 K only the OAGT values are given 
as recommended values. In iDR2 we disregarded the OAGT 
UVES values for Teff >6200 K. 


Note that, despite the large difference in resolution and 
spectral range, the comparison with benchmarks shows that 
the UVES Teff accuracy is only marginally better than GI¬ 
RAFFE’S, while log 5 and [Fe/H[ results from the two se¬ 
tups are of comparable accuracy. Our recommended values 
include Teff a nd [Fe/H[ for 11 stars and logg for 3 stars 
(see Sect. 5.7) from both the UVES and GIRAFFE setups. 
The comparison of our results for the same stars in the two 
setups shows that the Teff ratio (GIRAFFE/UVES) has a 
mean of 0.99 and a median of 1.00. The differences in [Fe/H[ 
(GIRAFFE-UVES) have a mean of 0.13 dex and a median 
of 0.16 dex. Among the 3 benchmark stars for which we 
give recommended \ogg from both GIRAFFE and UVES 
setups according to the criteria described in Sect. |5.7[ two 
are in the range of interest (log gw 4.0) and the maximum 
difference with the benchmark value is -0.09 dex. 
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Fig. 4. GIRAFFE benchmarks comparison. Black for [Fe/H]< 0.1, blue for [Fe/H]> 0.1. Size is inversely proportional to logg in 
the Teff plot, proportional to Ten in the logg and [Fe/H] plot. 

Table 3. Average differences from benchmark reference values (see Fig.and |^. The Arcetri’s [Fe/H] results are offset by 0.09 
dex before computing A [Fe/H] and CT[Fe/H[ (see text for details). 



(ATeff) 


(A log 5 ) 

(cr(logg)) 

(A [Fe/H]) 

(a([Fe/H[)) 

tact 


GIRAFFE 

OACT 

50. 

124. 

0.19 

0.29 

-0.04 

0.14 



OAPA 

18. 

120 . 

-0.15 

0.28 

-0.03 

0.16 



UVES 

OACT 

38. 

124. 

0.15 

0.26 

0.05 

0.18 



Arcetri 

55. 

95. 

0.14 

0.19 

0.10 

0.12 

0.00 

0.31 

CAUP 

34. 

96. 

-0.02 

0.28 

-0.03 

0.08 

0.04 

0.33 

UCM 

56. 

90. 

0.09 

0.25 

-0.01 

0.08 

0.04 

0.38 


Table 4. UVES results on the solar spectrum. 



ATeff 

trTeff 

Alogg 

cr log g 

A [Fe/H] 

cr [Fe/H] 

AC 


OACT 

-1. 

67. 

-0.14 

0.11 

0.06 

0.10 



Arcetri 

49. 

150. 

0.09 

0.20 

0.09 

0.07 

1.00 

0.15 

CAUP 

-48. 

59. 

-0.18 

0.12 

0.00 

0.07 

0.87 

0.08 

UCM 

9. 

48. 

0.00 

0.11 

0.03 

0.04 

0.75 

0.08 


5.7. Internal comparison 


The node-to-node comparison for the UVES individual 
clust er r esults before data selection and calibration (see 
Sect. 5.6) gives systematic differences in the ranges 80- 
160 K in. Teff, 0.1-0.3 dex in logg, and 0.06-0.17 dex in 
[Fe/H], while dispersions are in the ranges 160-260K in 
Tes, 0.1-0.3dex in logg, 0.13-0.45dex in [Fe/H]. The ap¬ 


plication of the data selection and calibration discussed 
in Sect. |5.6| reduces systematic differences below 100 K in 
Teff, and below 0.15dex in [Fe/H]. The final node-to-node 
mean dispersion in the recommended data is 110 K in Teff, 
0.21 dex in logg, and 0.10dex in [Fe/H]. These values are 
very close to the median dispersion: 106 K in Teff, 0.17 dex 
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Fig. 5. UVES benchmarks comparison. Black for [Fe/H]< 0.1, blue for [Fe/H]> 0.1. Size inversely proportional to logy in the T^s 


and ^ plot, proportional to T^s in the logy and [Fe/H] plot. Reference ^ from Bergemann et al. (2014, in prep.) (see also Smiljanic 
et al[^|20T4l. 


in logy, and 0.11 dex in [Fe/H], Biases in the recommended 
data are therefore successfully reduced. 

For the GIRAFFE results, systematic differences be¬ 
fore data selection are in the ranges 110-200 K in 0.4- 
0.8dex in logy, and 0.01-0.03 in [Fe/H], while dispersions 
are in the ranges 210-330K in Tgfj, 0.65-1.00dex in logy, 
and 0.17-0.26dex in [Fe/H], The situation here is more 
complex than in the UVES case. The problems to address 
are: 


(1) The OACT (RQTFIT) logy for PMS stars tend to be too 
high, clustering essentially on the 

(2) The OAPA logy for PMS stars tends to be too low, 
often lower than suggested by modelj^ 

(3) The PMS domain is contaminated by non-members with 
spurious logy in both logy-Teg diagrams. 

(4) The RGB in the OAGT logy-Teg diagra. m follows the 
calibrated relation taken from Gox (20001, while in the 
OAPA diagram it doesn’t. 

(5) The OAPA log y-T/ff diagram outside the MS, PMS, and 
RGB domains is sparsely populated, with both some 


® This is due to the basic criteria for defining the templates, 
identified as slow rotators, inactive stars and with no significant 
Li-absorption, which imply that no PMS star can be taken as 
template. 

An absolute calibration of the gravity-sensitive spectral in¬ 
dex in the PMS is very difficult (or impossible with currently 
available data) because of the lack of suitable PMS calibrators. 


very low and very high values, which are indicative of 
possible presence of some large errors. 

(6) The OAPA TeSi logy, and [Fe/H] are valid for usinf< 
90, 30, and 70 kms“^, respectively. 

(7) Because of a continuum normalisation problem on the 
Ha wings, in iDRl and iDR2 the OAGT parameters 
need to be discarded for Teff> 5500 K. 

In order to reduce biases as much as possible and pro¬ 
vide reliable recommended results we adopt the following 

solution: 

1. The OAPA are considered only for usinf < 90 
kms“^. The OAGT Tgg are considered only below 
5500 K. In cases where both the OAGT and OAPA Teff 
are available these are averaged. In all other cases the 
remaining value, if any, is adopted as recommended Teg. 

2. The OAPA [Fe/H] are considered only for vsini <70 
kms“^. The OAGT [Fe/H] are considered only below 
5500 K. In cases where both the OAGT and OAPA 
[Fe/H] are available these are averaged. In all other cases 
the remaining value, if any, is adopted as recommended 
[Fe/H], 

3. The OAPA logy are considered only for usinf < 
30 kms“^. The OAGT logy are considered only for 
Teg<5500K. In cases where both the OAGT and OAPA 
log y are available these are averaged if they differ by less 
than 0.3 dex. When they differ by more than 0.3 dex. 
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1.05 

1.00 

0.95 

^ 0.90 

0.85 

0.80 


Fig. 6. The OAPA gravity-sensitive 7 index vs. for all clus¬ 
ters analysed in iDRl and iDR2. The group of younger clusters 
(p Oph, Chal, and NGC2264) are clearly distinguishable from 
the group of older clusters (7 Vel and NGC2547). 


GIRAFFE PMS 



; O Cha I (2 Myr) 

-NGC2264 (3 Myr) 

; □ y Vel (5-10 Myr) 

- A NGC2547 (35 Myr) 


5000 4500 4000 3500 3000 

Teff (K) 


if the OACT logi; > 4.2 and the OAPA logg > 5, the 
OACT log <7 is given as recommended value. In all other 
cases we do not give recommended logg. 


4. The OAPA gravity-sensitive 7 index (Damiani et al. 


20141 is given as a recommended parameter for v sin i < 
MTms-h 


The application of such criteria leads to a final node- 
to-node mean dispersion in the recommended data of 98 K 
in Tqs, 0.23dex in log ( 7 , and 0.14dex in [Fe/H]. These val¬ 
ues are very close to the median dispersion: 95 K in Tgg, 
0.22dex in log 77 , and 0.14dex in [Fe/H]. Biases in the rec¬ 
ommended data are therefore successfully reduced in the 
GIRAFFE case too. 

When a recommended log 77 is not given, it may be 
still possible to identify an approximate evolutionary sta¬ 
tus based on the OACT and OAPA results. Those stars 
for which a trustworthy log 17 cannot be recommended are 
therefore flagged, when possible, as PMS, MS or post-MS 
stars according to the criteria listed in Table]^ 


Table 5. Criteria for the evolutionary status. 


Teff 

log ffOACT 

log 770APA 

Status 

< 5500K 

> 3 

3-4.2 

PMS 

< 5500K 

3-4.2 


PMS 



> 4.2 

MS 


< 3 

< 3 

post-MS 


The gravity-sensitive spectral index 7 obtained by the 


Damiani et al. (20141 approach can provide a rank order 
in age of the clusters. This can be seen, for the clusters 
analysed to date, in Fig.j^ where values for the younger 
clusters group (p Oph, Chal, and NGC2264) are clearly 
separated from those of older clusters group (7 Vel and 
NGC2547). However, both the scatter in 7 and the small 
age differences between clusters in the younger or the older 
group still prevent a clear separation in age. 
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5.8. Overview in the log g-Tgg plane 


As a final check on our recommended fundamental param¬ 
eters, we examine the log g-T^e diagram obtained with our 
data (Fig.[^ and compare it w i th the calibration of MK 


spectral classes reported in 

Gox (20001 and the theoretical 

PMS isochrones from 

Allarc 

et al. 

([2011 1 . We note the clus- 


tering of field stars on the MS and the RGB, as expected, 
while for the PMS clusters’ members a residual bias to¬ 
wards the MS and the RGB remains. In < 1/2 of the cases, 
log g values for PMS stars are located approximately where 
predicted by the models, although with large uncertainties. 



Fig. 7. logp-Tefi diagram for all targets. Grey filled circles are 
used for clusters’ non-members and stars not classified as GTTS 
nor WTTS. T he red dash ed lines are the dwarfs and giants se¬ 
quences from 


isochrones at 


Cox (2000 1 . The blue d ot-dashed lines are the 
., 5, 10, and 20 Myr from [Allard et ah (20111. 


5.9. Comparison between fundamental parameters derived 
from GIRAFFE and UVES 

A number of stars in the 7 Vel field have been observed 
with both UVES and GIRAFFE. For iDR2, our analysis 
produced and [Fe/H] values for 31 stars and log 77 for 16 
stars in this common sample. Note that the lower number of 
log g values is due to the application of the criteria described 
in Sect. |5.7| which were applied to iDR2 but not to iDRl. 
The comparison of the recommended values for this sample 
is satisfactory (see Fig. 1^ for iDR2) and support the valid¬ 
ity of our approach both in the parameters determination 
and in the derivation of the recommended values. A similar 


comparison is reported in Spina et al. (2014bI for iDRl. 
Indeed the reproducibility of the parameters obtained with 
the higher resolution and larger wavelength coverage from 
UVES using a much smaller wavelength range and a lower 
resolution as in GIRAFFE is a remarkable achievement and 
increases our confidence in our parameters determination 
from the much larger GIRAFFE sample. 


5.10. Veiling vs. Ha emission 

For iDRl, the RQTFIT veiling parameter was adopted as 
the recommended one. In iDR2, however, it was recog¬ 
nised that some residual nebular emission remained after 
sky-subtraction, particularly in NGG2264, which were not 
sufficiently masked in the RQTFIT calculations. As a con- 
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UVES UVES UVES 

Fig. 8 . Comparison between fundamental parameters derived from GIRAFFE and UVES spectra in the 7 Vel held (iDR2). 


sequence, the RQTFIT veiling parameter for NGC2264 was 
clearly overestimated and the OAPA solution was adopted 
as recommended in iDR2. Note that this does not invali¬ 
date the results of iDRl as 7 Vel and Chal spectra are not 
affected by residual sky emission in the reduced spectra. 


Frasca et al. (20141 found a positive correlation between 
Ha ffux and r in the iDRl data for Chal objects with 
r > 0.25, for which the Spearman’s rank analysis yielded 
a coefficient p = 0.58 with a significance of cr = 0.003. 
The same analysis for all clusters in iDR2 gives a coeffi¬ 
cient p = 0.39 with a significance of u = 0.004. However, 
a correlation between r and Ha 10% or M^(Ha) is not evi¬ 
dent in the iDR2 data (see Fig.l^, where we do see an in¬ 
crease of the upper envelope witheither Ha 10% or W (Ha), 
but the large scatter makes the correlation not significant. 
This is at variance with w hat expected from previous work 
(e.g.. White fc Basri|20'03 1 and therefore it outlines possible 
limitations in our veiling determination. Further validation 
based on comparison with different methods is deferred to 
future work. 


6. Derived parameters 

6.1. Li abundance 



recommended log e(Li) is derived using the recommended 
fundamental parameters and recommended W (Li) as input. 
Uncertainties were obtained by propagating the Teff and 
lU(Li) uncertainties. 

The approach adopted in the GIRAFFE case has the 
advantage of allowing us to focus on the accuracy of the fun¬ 
damental parameters and IF (Li), relying then on the best 
GOG available to derive node-specific and recommended 
loge(Li). Note that the two GOGs adopted do not join 
smoothly at 4000 K, but the interpolation scheme ensures 
a smooth transition between the two regimes. A derivation 



200 300 400 500 600 700 

Ha 10% 



10 


100 

W(Ha) 


Fig. 9. Veiling parameter r vs. Ha 10% (top panel) and vs. 
IF(Ha) (bottom panel) for iDR2. 


of a self-consistent GOG in the whole range is planned 
as a future improvement. T he n ode-to-node dispersion in 
the GIRAFFE case (see Fig.jToj for the whole iDR2) then 
propagates only from the lF(Ti) measurements and shows 
a fairly random distribution wit a median of O.lTdex. 
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In the UVES analysis, the OACT and Arcetri loge(Li) 
were derived as in the GIRAFFE case. The CAUP and 
UCM nodes, on the other hand, derived loge(Li) by a stan¬ 
dard LTE analysis using the driver abfind in the revised 
version of the spectra l sy nthesis code MQQG (Sneden 19731 
(see also Sects. and . CAUP used the 2010 version of 
MQQG, while UCM used the 2002 and 2013 versions for iDRl 
and iDR2, respectively. Uncertainties were estimated by 
varying each atmospheric parameter within its uncertainty 
range to derive the propagated uncertainty in log e(Li). The 
propagated uncertainties where then combined quadrati- 
cally. In this case the recommended value is given as the 
average of all nodes’ estimates available with a cr-clipping 
when at least three measurements are available. FigurefTO] 
shows that also in this case the node-to-node dispersion 
have a fairly random distribution, with a median uncer¬ 
tainty of 0.12 dex. 

Possible ®Li contribution was neglected in all cases. 


Li I UVES 



Li I GIRAFFE 



6.2. Other elemental abundances 

Elemental abundances were computed by three nodes 
(Arcetri, CAUP, and UCM) when good quality UVES spec¬ 
tra were available in stars not affected by veiling and/or 
large usini. 

The Arcetri no de computed abunda nces using FAMA. We 
refer the reader to Magrini et al. (2013) for a description of 
the method and the way in which lines are selected for the 
abundance analysis. 

The CAUP node derived individual abundances using 
the driver abfind in the 2010 version of MQQG (see [Neves 
|et al. 2009| Adibekyan et al. 2012 for details) and equiv¬ 
alent wTHthslneasured^witntlieARES code. The line list 
for elements other than Fe (with atomic number A <28) 
was selected through the cross-matc hing between the line 
list used by Adibekyan et al. (2012|) and the line list pro¬ 
vided by Gaia-ESO. F'or elements with A >28, lines that 
were suitable for W measurements (as tested by the Gaia- 
ESO line-list working group) were first selected and from 
these the ones that ARES was able to measure were used. 
The atomic data from the Gaia-ESO Survey was adopted. 
CAUP considered hyperfine splitting in the analysis of Cu, 
Ba, Nd, Sm and Eu abundances, i.e. for all the elements af¬ 
fected with A > 28 (using the driver blends in MQQG). The 
errors of the abundances is given as the line-to-line scatter 
(when more than one line is measured). 

The UCM node adopted an approach similar to CAUP. 
For iDRl, two line-lists were prepared: one for dwarfs 
(logg > 4.0) and one for giants (log 5 < 4.0). For iDR2 
five line lists were used as done for the stellar parameters 
(see Sect.|^. A total of 13 elements were analysed: Fe, the 
a-elements (Mg, Si, Ca, and Ti), the Fe-peak elements (Cr, 
Mn, Co, and Ni), and the odd-Z elements (Na, Al, Sc, and 
V). To obtain individual abundances, the equivalent widths 
are fed into MQQG and then a 3(T-clipping for each chemical 
element was applied. 

The elements for which at least two nodes derived abun¬ 
dances for at least one star and that were considered in the 
recommended results are: Na, Mg, Al, Si, Ca, Sc, Ti, V, 
Cr, Mn, Fe, Co, Ni, Zn, Zr, Mo, Ce. Only one node re¬ 
sults were considered for: Cu, Y, Ba, La, Pr, Nd, Sm, Eu. 
Abundances are from the neutral species except for Ba, La, 
Ce, Pr, Nd, Sm, Eu, for which they are from the ionised 
species. 

The node-to-node dispersions of elemental abundances 
in iDR2 is shown in Fig. 11 In general, « 90% of the results 
for each elements have dispersions below « 0.2 dex. How¬ 
ever the tail of the distributions extends to higher values in 
more difficult cases for which differences that arise from the 


Fig. 10. Li abundance uncertainty histogram for all sources 
in iDR2. For the GIRAFFE spectra the loge(Li) uncertainty 
is propagated from the uncertainty in Tes and IT (Li). For the 
UVES spectra the node-to-node dispersion is considered. A solid 
line is used for the cumulative probability (right ordinate axis). 
See text for details. 


different W measurements (see Smiljanic et al. 2014) and 
line selection strategies play a role. The dispersion tends to 
be higher also for abundances of ions like Tin and Cm. 
Poor agreement is found for Zni and Zrll. Note that abun¬ 
dances for elements which require hyper-fine splitting were 
provided by the CAUP node only. 

The internal precision is comparable with that of the 
UVES spectra of FGK-typ e analysis (excluding stars in the 
field of young open clu ster; Smiljanic et al. |2014[ ) . Note that 


Smiljanic et al. (2014) make use of the median of the ab- 
solute deviations from the median of the data (MAD) to 
quantify the node-to-node dispersion, but this cannot be 
used here because of the small number of nodes providing 
abundances. The dispersion from the mean used here should 
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overestimate the node-to-node dispersion with respect to 
the MAD, although this is mitigated by the cr-clipping ap¬ 
plied. Overall, all this indicates that our internal precision 
for elemental abundan ces is roughly at the same level of the 

Smiljanic et al. (20141 one. _ 

A survey inter-comparison with Smiljanic et al. (20141 


et al. (2014). They found discrepancies of ^ 0.8 dex for Cha 


1 and ~ 0.7 dex for 7 Vel on average. Frasca et al. (2014) 


results on the common calibration open cluster NGC6705 
was carried out for all elemen ts except Ce, La, Pr, Sm, for 
which results did not pass the Smiljanic et al. (20141 quality 
control criteria. The inter-comparison was satisfactory and 
confirm ed the comparable precision with the [Smiljanic et ^ 
(2014) results. 


also compared the results obtained for Chal with literature 
values, finding a fair agreement, with differences that can 
be ascribed to variability, different methodologies and the 
use of different evolutionary models. 

Mass accretion rates derived for all clusters in the first 
18 months of observations vs. stellar mass are shown in 
Fig.IT^ Stellar mass is estimated from the recommended 


Teff and the age of the cluster using the Baraffe et al. (1998 1 


63. Mass accretion rate 


-6 


(e.g., Muzerolle et al. 

2005 

Herezeg & Hillenbrand 

2008 

Alcala et al. 2014|). As for t 

he 7 Vel and Chal cases dis- 

cussed in Frasca et al. 

(2014), however, the large scatter in 
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M prevents us to make a meaningful comparison with such 
a relationship. The Spearman’s rank correlation analysis 
for Chal gives, in the iDR2 case, p — 0.43 and u = 0.005 , 
i.e. a higher significance than found by Frasca et al. (2014) 


- 0.4 


- 0.2 


0.0 

Log M 


0.2 


0.4 


Fig. 12. Mass accretion rates vs. mass for all clusters in iDR2. 
Symbols and colours as in Fig.j^ The dashed line represents the 
M oc relationship. 


in the iDRl case (p = 0.26, cr = 0.16 for M derived from 
Ha 10%), which indicates a better accuracy of our recom¬ 
mended iDR2 Ha 10% parameter. Amongst the younger 
clusters in our sample, we find p = 0.47 and a — 0.14 for 
p Oph, while the correlation is rather poor for NGC2264 
(p = 0.19, a = 0.022) possibly because of the larger uncer¬ 
tainties due to the residual nebular emission in the spectra 
of this cluster. Note that the scatter in Fig.[^is dominated 
by NGC2264. Ignori ng this cluste r , the s catter is consistent 
with what found b y Alca la et al. ([2014 ) in their validation 
of the Natta et al. ( 2004p relationsnip. Interestingly, for the 
older clusters in our sample we find a not significant cor¬ 
relation in 7 Vel (p = 0.29, a = 0.247) but a well defined 
correlation in NGG2547 (p = 0.89, a = 0.018). In both 
such cases, two kinematically dist inct populations with dif¬ 
ferent ages have been discovered (Jeffries et al.|20f4 Sacco 


et al. 


2015), whose possible consequences in the M vs. M 


relationship still need to be explored. 


Mass accretion rates are estimated from the Ha 10% Chromospheric Ha and H/3 flux 


using the Natta et al. (2004 Eq. (1)) formula. The use of 
alternativ e methods, i.e mak ing use of the IF(Ha), is dis¬ 
cussed in 
the Gaia- 


Frasca et al. (2014) and will be implemented in 
USO P MS analysis in futur e data releases. 

The use of the Natta et al. (2004) relationship has, un¬ 
doubtedly, the advantage of allowing a simple estimate of 
Mace from just the Ha 10%. The accuracy and validity of 
this empi rical relationship has , however, been questioned 
(see, e.g., Gostigan et al. [2012 and references therein), es¬ 
pecially in cases when only single epoch observations are 
available. 


After the ROTFIT determination of the fundamental pa¬ 
rameters, a best matching template within the library of 
slowly-rotating inactive stars is identified. The chromo¬ 
spheric excesses AkF(Ha)chr and A1F(H/J)c hr are derived 
using a spectral subtraction method (see, e.g., Bardenjl985 


Recently, Alcala et al. (2014) have computed the accre¬ 
tion rate by modelling the excess emission from the UV to 
the near-IR and provided empirical relationships between 
accretion luminosity and the luminosity of 39 emission lines 
from X-Shooter spectra. In particular, they have shown that 
the comparison between Mace derived through primary di- 
agnostics (like the UV-excess) and that obtained with the 
Natta et al. (2004) relationship has a large scatter, with 
this latter tending to underestimate Mace for Ha 10%< 400 
kms“^ and to overestimate Mace for Ha 10%> 400 kms“^. 

A comparison of the iDRl Mace of 7 Vel and Ghal with 
mass accretion ra tes derived from line luminosity and the 
Hartmann (1998) relationship has been presented in Frasca 


[Frasca fc: Gatalano|1994| [Montes et al.|1995| and references 
therein) that has been extensively used in the past. The 
photospheric flux is removed by subtraction of the spectrum 
of an inactive template star with very close fundamental pa¬ 
rameters, rotationally broadened at the target usini, over 
the line wavelength range. Such chromospheric W excesses, 
AlF(Ha)chr and AkF(H/3)chrj are then converted to flux, 
F(Ha)chr and F(H/3)chr, by multiplying it by the t heoreti¬ 
cal continuu m flux at the line’s wavelength (see, e.g., Frasca 
et al.[[20T4 and references therein). It may be argued that 


even the templates may have some chromospheric basal flux 
(see, e.g.. Judge fc Garpenterjl998[ and references therein), 
also variable in time following the stellar cycles (see, e.g., 
Schroder et al. 12012 ) which a de tailed semi-empirical NLTE 


chromospheric m odelling (e.g., Houdebine & Panagi 1990 
Lanzafame 1995) could take into account. 'This latter is. 


The results shown here are for a mixing length parameter 
a = 1.5; in this analysis, however, the choice of a is uninfluential. 
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however, unpractical for applications to large datasets like 
the Gaia-ESO one. Furthermore, the chromospheric flux in 
young stars is much larger than the basal flux, so that this 
latter can be safely neglected. 

Results for 7 Vel and Chal (iDRl) are discussed in 


Frasca et al. 


( 2014[ ), who were able to discriminate be- 
tween chromospheric-dominated and accretion-dominated 
Ha flux. AbF(Ha)chr vs. T^s for all clusters observed in the 


first 18 months of observations (iDR2) is shown in Fig. 13 
We note that the chromospheric activity - accretion on 


ual nebular emissio n. This dividing line was also found by 
Frasca et al. (2014|) to be in remarkable agreement with the 
saturation lirnit adopted by jBarrado y Navascues & Martini 
(|2^ to separate CTTS and WTTS. 


viding line proposed by Frasca et al. (2014) (log Fna = 
6.35 -I- 0.00049(reff — 3000)) delimits quite neatly the two 
regimes in this larger sample as well, with some larger 
uncertainties in the case of NGC2264 likely due to resid¬ 


Finally, in Fig. 14 we show AkF(Ha)chr vs. usini for 
all young clusters observed in the first 18 months of ob¬ 
servations. While a full discussion on the activity-rotation 
relationship is deferred to future work, we note that our 
data display a Tgff-dependent activity-rotation correlation 
regime at low vsini, followed by a Teff-dependent satura¬ 
tion regime at high usini, as expected. The behaviour at 
different is quite neatly distinguishable, which further 
confirm the overall consistency of our results. 


Article number, page 16 of 


22 


















































































A. C. Lanzafame et al.: Gaia-ESO Survey: The analysis of pre-main sequence stellar spectra. 



Fig. 13. Chromospheric Ha flux vs. for all young clusters 
observed in the first 18 months of observations. Symbols and 
colours as in Fig.lH with filled (open) symbols used for CTTS 
(WTTS). The dashed line represents the chrom ospheric activity 
- accretion dividing line of Frasca et al. (|2014 1. 
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Fig. 14. Chromospheric Ha equivalent width excess vs. vsini 
for all young clusters observed in the first 18 months of obser¬ 
vations. Colour coding is used for T^g. 


7. Summary and conclusions 

The Gaia-ESO PMS spectrum analysis provides an exten¬ 
sive list of stellar parameters from spectra acquired in the 
FLAMES/GIRAFFE/HR15N and FLAMES/UVES/580 
setups in the field of young open clusters. These include raw 
parameters that are directly measured on the input spec¬ 
tra (VF(Ha), Ha 10%, and W (Li)), fundamental parameters 
(Teff) logg, [Fe/H], usinf, and r), and derived parame¬ 
ters (loge(Li), loge(X)), Mace, AVT(Ha)chr, AVF(H/3)chr, 
F(Ha)chr, and F(H/3)chr) which require prior knowledge of 
the former. Our analysis strategy is devised to deal with 
peculiarities of PMS stars and young stars in general such 
as veiling, large broadening due to fast-rotation, emission 
lines due to accretion and/or chromospheric activity, and 
molecular bands. The analysis is also made robust against 
residual sky-background or foreground features that can¬ 
not be completely removed as in the case of inhomogeneous 
nebular emission. 

The availability of different methods for deriving stel¬ 
lar parameters increases the confidence on the output of 


our analysis. It allows us to efficiently identify and discard 
outliers, like those deriving from failed fits or problems in 
the input spectra, as well as deriving realistic uncertainties 
from the internal dispersion of the data. For Teff and log g 
the external precision is estimated by comparison with re¬ 
sults from interferometric angular diameter measurements. 
These are estimated to be « 120 K r.m.s. in and «0.3 
dex r.m.s. in logt/ for both the UVES and GIRAFFE se¬ 
tups. The comparison with Teff derived from photometry 
for a selected group of stars in 7 Vel with the same fore¬ 
ground extinction and free from accretion signatures gives 
an agreement of « 260 K r.m.s. Our recommended [Fe/H] 
results agree with assessed literature values for such a set 
of benchmark stars within «0.15 dex r.m.s. A comparison 
with previous [Fe/H ] determination for Gha I is discussed in 
Spina et al. (2014a I. Weakness or limitations of the meth- 
ods used were identified by the node-to-node comparisons 
and by comparison with benchmark stars. 

The observation strategy poses significant challenges to 
the analysis, since, for optimising the observation time, 
most of the relevant observations are carried out in just the 
FLAMES/GIRAFFE/HR15N setup. For our purposes, the 
wavelength range of this setup is the best available in the 
optical, as it contains very important diagnostics for young 
stars like the Ha and Li line. At the same time, surface 
gravity diagnostics in the HR15N setup are poorer than in 
other wavelength ranges and still not modelled with suffi¬ 
cient accuracy. Teff determination for spectral types earlier 
than early-G is also challenging since it is based mostly on 
the Ha wings. For such a wavelength range, two methods 
based on the comparison with spectra or spectral indices of 
template stars have proved effective in providing fundamen¬ 
tal parameters. A satisfactory self-consistency of the results 
have been achieved, at the expense of discarding logg val¬ 
ues when a sufficient agreement between the two methods 
cannot be reached. In such cases, however, it is still possi¬ 
ble to provide an evolutionary flag, as it can be established 
with confidence whether the star is in a PMS, a MS or a 
post-MS stage. An uncalibrated gravity-sensitive spectral 
index is also provided, useful for a rank order in age. 

The reproducibility of the parameters obtained with 
the higher resolution and larger wavelength coverage from 
UVES using a much smaller wavelength range and a lower 
resolution as in the GIRAFFE/HR15N setup, together with 
the comparable accuracy and precision achieved in the two 
setups, is a remarkable achievement of this work. This al¬ 
lows us to provide with confidence parameters for the much 
larger GIRAFFE sample. 

The Gaia-ESO is an ongoing project and this paper 
describes the PMS spectrum analysis carried out on the 
first two data releases. Work is ongoing to improve fur¬ 
ther our analysis for the next releases. The tables with 
the public release results will be available through the 
ESO data archiv^^ and through the Gaia-ESO Survey sci¬ 
ence archiv^^ hosted by the Wide Field Astronomy Unit 
(WFAU) of the Institute for Astronomy, Royal Observatory, 
Edinburgh, UK 
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Fig. A.l. Parameters of the whole set of ROTFIT templates 
adopted for the Gaia-ESO analysis. 


Appendix A: ROTFIT templates 


The method implemented in ROTFIT relies on existing de¬ 
termination of fundamental parameters for the template 
stars and, in the past, the PASTEL catalogue (Soubiran 


et al. 20101 has been used as input. To ensure homogene¬ 


ity amongst the Gaia-ESO spectrum analysis, however, the 
templates’ parameters have been re- determined using Fast 
Automatic MOOG analysis (FAMA, Magrini et al. 20131 
adopting the Gaia-ESO recommended model atmospheres 
and atomic parameters. The templates’ parameters were 
updated for most of the stars in the range from mid-F 
to late K, while for the M dwarf we adopted the pa ram¬ 
eters recently d e termine d by Rojas-Ayala et al. (2012) and 
Boyajian et al.| (|2012b|). The parameters of the whole set 
of ROTFIT templates adopted for the Gaia-ESO analysis is 
illustrated in Fig. EH In Fig. |A.2| we show the compari¬ 
son between the PASTEL and FAMA parameters. Average 
differences are 30K, O.lOdex, and 0.04dex for Tes, logg, 
and [Fe/H], respectively. Standard deviations are 123 K, 
0.27dex, and 0.11 dex for Tgg, logg, and [Fe/H], respec¬ 
tively. A table with the ROTFIT templ ate parameters used 
for the Gaia-ESO analysis is reported in Frasca et al. (20141. 


Appendix B: ROTFIT masks 

Spectra of accreting stars or embedded in a dense cloud 
requires wavelength masks to exclude the residual nebular 
emission features that still remain after the data reduc¬ 
tion process. The Ha profile in non-accreting young stars 
is nonetheless affected by a significant chromospheric con¬ 
tribution, which also varies in time. The Ha core must 
therefore be masked out to avoid considering the non- 
photospheric contribution to the line profile. On the other 
hand, the Ha wings are essential fundamental parameters’ 
diagnostics, especially in GIRAFFE/HR15N spectra of F- 
and G-type stars for which the rest of the pass-band offers 
very poor constraints. In some cases (e.g., old and inac¬ 
tive stars in the field) the whole of the Ha profile can be 
used. Because of the Li depletion occurring in the stellar 
interior, W (Li) decreases rapidly with age in PMS stars of 
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Elodie Templates: FAMA vs PASTEL 




Elodie Templates: FAMA vs PASTEL 



Fig. A.2. Comparison between PASTEL and FAMA fundamental 
parameters of the ELODIE templates used by ROTFIT. See text 
for details. 


later spectral type and can be very different in stars with 
similar fundamental parameters; therefore this line must 
also be masked out. The Ha and Li (6707.84 A) lines must 
be masked in the spectra of accreting stars too, but in such 
cases the Ha mask must be wide enough to include the 
wings of the lines that can be very broad. 

Therefore, the measurement of the raw parameters as 
a first step in the analysis process allows us to divide the 
spectra into three classes: 

— NHL. Spectra with negligible Li absorption and no Ha 
emission for which only a narrow Ha-core mask (±2 A) 
is required and no r evaluation is carried out; 

— HL. Spectra with significant Li absorption, Ha core 
emission, or both, for which accretion, and therefore r 
evaluation, can be excluded, but require a slightly larger 
Ha-core mask (±5 A); 

— HLV. Spectra with accretion signatures, which require 
a mask for the entire Ha profile (±20 A), plus the eval¬ 
uation of r. 

In the case of HL and HLV classes, a mask of ±3 A is applied 
around the Li line core. 


Appendix C: Comparison with from 
photometry 

A more extensive comparison can be made, at this stage, 
with Teff derived from photometry for cluster’s members 
that have all the same foreground extinction and that are 
not significantly affected by colour excesses due to the pres¬ 
ence of circumstellar material. Fortunately, this turned out 
to be feasible for 7 Vel, as the extinction is fairly uniform 
and many likely members are free from large colour excess. 
BVI photo metry of 7 Vel was p resented in |Pozzo et al.| 


(20001 and in Jeffries et al. (2009). The photometry was 
taken in the Harris H, V and Kron-Cousins I filters and was 
converted to the standa r d Joh nsons-Cousins photometric 
system by Jeffries et al. (20( ^. This dat a have been sup- 
plemented with 2MASS (Cutri et al. 20031 and Spitzer data 
( Hernandez et al(]|2008 l where available. 

We used two different methods for deriving Teff from 
photometry, one based on the {V — I) vs. {B — V) colour- 
colour diagram (Tqc), the other one based on the simul¬ 
taneous fit of all available magnitudes, from the optical to 
the Spitzer bands (Tsed)- 

The iy — I) vs. (B — V) colour-colour diagram for likely 
members of 7 Vel, corrected for th e foreground ext inctio n as 
estimated by Jeffries et al. (20091, is shown in Fig. C.l The 
membership is based on Li and radial velocity and about 


85% of the stars i n Fig. C.l are expect ed to be actual clus¬ 
ter members (see Jeffries et al. 2014). The narrowness of 
this locus indicates that there is little differential extinc¬ 
tion or veiling in 7 Vel. When excluding stars with radial 
velocities larger than 9 kms“^ from the mean velocity of 7 
Vel, this locus narrows even further. As Fig. |C.l| shows, nei- 
ther the empirical Teff-to-co lour conversion for the ZAMS 
(Kenyon & Hartmann 1995 hereafter KH95) nor th e theo¬ 
retical HT- Settl ZAMS and 7 Myr PMS isochrones (Allard 
et al.|[ 20 TT l overlap with the observed locus in the colour- 
colour diagram. 

The theoretically computed colours from BT-Settl do 
not reproduce accurately the stellar magnitudes possibly 
because of a still incomplete description of the opacity. 
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Fig. C.l. V — I versus B — V colour-colour diagram for likely 
members of 7 Vel. Stars with radial velocities larger than 9 
kms“^ from the mean velocity of 7 Vel are crossed out. The em- 
pirical main sequence locus relationship from [Kenyon fc Hart^ 
(19951 is over-plotted in green, the theoretical main se¬ 


quence locus from BT-Settl (Allard et al. 20111 in cyan, the 
theoretical pre-main sequence locus for an age of 7 Myr from 
BT-Settl in purple, and finally the locus obtained from Eq. (C.l I 
in red. All the colour loci plotted as lines h ave been shifted as- 
suming an extinction ot E{B — V) = 0.038 ([Jeffr ies et a O2009 ( 
and E{V — I)/E{B — V) — 1.6 ( Rieke fc Lebofs^y|T^^^P 


However, we can use them to define a new locus based on 
the assumption that the models correctly predict the colour 
difference due to gravity effects. If this assumption holds, 
we can add the colour shift between the theoretical ZAMS 
and PMS isochrones to the empirical ZAMS colours to get 
the actual PMS intrin sic colours, a pro cedure conceptually 
similar to that used by 
temperature, a new co. 


Bell et al. (2013). For every effective 
our has been calculated using: 


(B-V) 


{B — P)kH 95.ZAMS + (C.l) 

{B — P)bt-s,pms — {B — P)bt-s,zams 


with a similar equation for {V — I). This new locus clearly 
describes the observed locus in 7 Vel much better than 
the alternatives (see Fig. C.l[ ). The fit can probably be im¬ 
proved further by varying the age an d extinction within the 
ranges given by Jeffries et al. (20091, but we do not explore 
that here. 

Assuming a con stant extinction of E{B — V) = 0.038 
( Jeffries et al.|[2009 l and zero veiling, which is likely to be 
approximately true for most stars lying on the observed nar¬ 
row locus, we can fit Tcc to the observed colours through 
a least-square minimisation, which has been generalised for 


two dimensions with correlated uncertainties (Hogg et al. 


2010 ): 


Xcc = ( ^BV Ap/ ) 


0-y CTy + erf 


Asy \ 
Av/ ) 

(C.2) 


where /S.bv = A(i3 — V) and Ayi = A(V — I) give the 
difference between the observed colours and those predicted 
from our new locus as a function of the temperature of 
the star. These colour differences are multiplied with the 
inverse of the covariance matrix, which is given here as a 


function of the photometric uncertainties (Tb, uy, and cr/. 
The temperatures computed in this way are based on the 
KH95 colour-temperature conversion of ZAMS stars with 
the colours adjusted for the lower surface gravity of PMS 
stars. 


In Fig. C.2 (top right panel) we compare the recom¬ 


mended Teff with Tec for likely 7 Vel members. A further 
selection has been applied by considering Tec with Xcc 
< 7 (i.e., consistent with being drawn by chance at the 99% 
level) to avoid considering stars that may be significantly af¬ 
fected by colour excess due to the presence of circumstellar 
material. Note that the formal uncertainties of Tec turned 
out to be excessively small as a result of the definition of 
the locus in the (V — I,B — V) plane and the insensitiv¬ 
ity of Tcc to B — V at low temperatures. The comparison 
shows that the agreement is mostly within the estimated 
uncertainties, although the spectroscopic Teff is systemat¬ 
ically higher than Tcc above « 3600 K, with an average 
difference of « 180 K and standard deviation a « 250 K. 

As an alternative approach, Teff from photometry has 
also been derived by taking all photometry available from 
optical, 2MASS, and Spitzer into account. In this case we fit 
Fsed by a- downhill sim plex multidimensional minimisation 
( Nelder fc Mead||l965 l of 


XSED 


= E 


Xi — Wi 
<Ji 


(C.3) 


where Xi = (V — M\) is the observed colour derived 
from the V-band magnitude and each of the photomet¬ 
ric magnitudes available {M\), (Ji is its uncertainty, and 
Wi = Wf(Tefr, logff , [Fe/H]) is the theo retical colour from the 
BT-Settl models (Allard et al. 2011) that are interpolated 
in Teff and log 5 with |Fe/H| fixed at the solar value. Stan¬ 
dard deviation has been estimated using Monte-Carlo sim¬ 
ulations with 1000 random synthetic realisations for each 
star. This method is, effectively, a fit of the stellar spectral 
energy distribution (SED) and has the advantage of con¬ 
sidering colours more sensitive to Teft than just {V — I) and 
{B — V) in the temperature range of interest. The weakness 
of this method lies mostly in the theoretical model that, al¬ 
though being amongst the most advanced available to-date, 
yet do not accurately reproduce observed PMS colour (see. 


e.g. 


Bell et aO2012 Stauffer et al.|[2007 ). 


Figure C.2 (top left panel) shows the comparison be¬ 


tween Tsed and Tcc for all GIRAFFE’S spectra in the 7 
Vel field for which < 7 (no significant IR-excess), ir¬ 
respective of Li- and RV-membership. Although some sys¬ 
tematic deviations in some temperature ranges are present 
(viz., Tcc <3800K), mainly due to the way in which Tcc 
is derived, the two models agree within the error bars, the 
mean difference is ~ —40K, and the standard deviation 
(T ~ 150 K. The comparison of the recommended Teff and 
Tsed for GIRAFFE is shown in the bottom left panel of 
Fig . |C.2| The two sets generally agree within the error bars, 
the mean difference is ~ 180 K and a ~ 240 K. Finally, the 
comparison of the UVES recommen ded Teff with Tsed for 
7 Vel radial velocity members (as in Spina et al.||2014b | is 
shown in the bottom right panel of Fig. [C.2| 

Note, finally, that the intrinsic variability of the tar¬ 
gets together the non-simultaneity of the spectroscopic and 
photometric observations must also play a role in the com¬ 
parisons presented here. Considering also uncertainties and 
likely spread in age, which are not reliably estimated as yet. 
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Fig. C.2. Top left panel: comparison of Tsed and Tec for all GIRAFFE’S spectra in the 7 Vel field for which Xcc < irrespective 
of Li- and RV-membership. Top right panel: recommended Teff vs. Tec for GIRAFFE’S spectra of likely members of 7 Vel. Bottom 
left panel: recommended Teff vs. Tsed for GIRAFFE’S spectra of likely members of 7 Vel. Bottom right panel: recommended Tefj 
vs. Tsed for all UVES’s spectra in the 7 Vel field. 


the comparison with photometry can be considered quite 
satisfactory. 
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